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TABLES 


4340  Steel  Contained  Fragmenting  Cylinder  Experimental 
Parameters  and  Qualitative  Damage  Results  . 

Material  Properties  for  Several  Steels  . 


1  SUMMARY 


The  purpose  of  this  program  was  to  improve  the  previously  developed 
computational  model  for  adiabatic  Shear  Band  Nucleation  And  Growth  (SNAG) . 
Specifically,  the  goals  were  to  perform  careful  experiments  to  obtain 
quantitative  data  for  improving  the  growth  law  for  the  shear  bands,  to 
derive  and  test  a  theory  for  the  threshold  conditions  for  nucleation  of 
shear  bands,  and  to  establish  optimal  experimental  procedures  for  measuring 
the  material  properties  governing  shear  band  behavior. 

In  pursuit  of  these  goals,  we  performed  seventeen  contained  frag¬ 
menting  cylinder  experiments  with  three  hardnesses  of  4340  steel.  (See 
Figure  1  and  Table  1).  Quantitative  data  for  the  number,  size,  and 
orientation  distributions  of  shear  bands  for  the  Rc40  experiments  were 
obtained  (See  Figures  5  through  12);  the  harder  materials  were  found  to 
be  more  susceptible  to  shear  banding.  Preliminary  correlations  of  these 
data  with  computer  computations  of  the  strain  and  stress  histories  yielded 
first  estimates  of  the  growth  parameters  in  the  SNAG  model  (Figures  15 
through  19).  A  strong  dependence  on  strain  rate  is  apparent. 

A  theory  for  the  threshold  nucleation  conditions  for  shear  bands 
was  derived  and  tested  against  experimental  data  (Figure  32).  The 
theory  predicts  that  the  shape  of  the  work  hardening  curve  is  the  key 
material  property  governing  resistance  to  shear  banding.  The  comparison 
with  experimental  data  suggests  that  the  postulated  threshold  criterion 
is  promising,  but  more  comparisons  are  needed  to  confirm  its  validity. 

Several  experimental  test  procedures  for  measuring  shear  band 
behavior  were  studied.  It  was  concluded  that  the  contained  fragmenting 
cylinder  test  used  in  this  program  is  the  best  single  test  procedure 
for  measuring  the  SNAG  parameters.  However,  for  specific  applications 
such  as  armor  penetration,  supplementary  tests  that  approximate  the 
geometries  of  interest  are  desirable. 
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The  results  of  the  program  confirm  that  future  work  should  continue 
the  quantitative  data  analysis  of  the  contained  fragmenting  cylinder 
experiments,  as  well  as  further  computational  analysis  to  more  clearly 
define  the  dependence  of  the  SNAG  parameters  on  key  material  properties. 
The  resulting  understanding  of  the  shear  banding  process  in  A3A0  steel 
will  allow  extrapolation  to  other  materials  such  as  those  used  for  armor 
and  penetrators. 
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2  BACKGROUND  AND  APPROACH 


Shear  banding  refers  to  the  development  of  internal  regions  of 
localized,  intense  plastic  strain.  It  is  the  governing  failure  mechanism 
in  many  cases  of  armor  penetration  by  plugging,  in  the  disintegration 

of  long  rod  penetrators,  and  in  the  fragmentation  of  exploding  rounds. 

12 

In  previous  work  for  the  Naval  Surface  Weapons  Center,  the  Ballistics 
Research  Laboratory,  and  the  Defense  Advanced  Research  Projects  Agency, 
SRI  International  performed  experimental  and  analytical  work  to  develop 
a  computational  model  to  describe  and  predict  in  quantitative  detail  the 
shear  banding  process.  The  resulting  model,  called  SNAG  for  fUiear 
Nucleation  And  Growth,  was  used  to  make  reasonable  predictions  of  the 
fragmentation  behavior  of  cylinders  of  Armco  iron,  HF-1  steel,  and  an 
actual  HF-1  steel  fragmenting  projectile. 

The  overall  objectives  of  the  present  program  have  been  improve  the 
SNAG  model  for  nucleation,  growth,  and  coalescence  of  adiabatic  shear 
bands,  and  to  establish  optimal  procedures  for  characterizing  the  shear 
band  behavior  of  materials.  Our  previous  work  showed  that  the  nucleation 
rate  and  coalescence  process  are  fairly  well  modeled  once  the  bands  have 
begun  to  form;  the  areas  of  most  uncertainty  are  the  threshold  conditions 
that  trigger  nucleation,  and  the  equations  that  describe  the  shear  band 
growth  velocities.  The  specific  objectives  of  the  present  program  were 
therefore  to: 

1.  Obtain  quantitative  experimental  data  to  improve  the 
nucleation  and  growth  equations 

2.  Derive  the  threshold  conditions  for  shear  band  nucleation 

3.  Establish  opt imal  experimental  procedures  for  measuring  the 
material  properties  governing  shear  band  behavior. 
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In  the  design  of  fragmenting  rounds,  the  growth  and  coalescence 
equations  are  especially  important  because  they  determine  the  fragment 
size  distributions.  Hence,  our  first  specific  objective  is  to  obtain 
quantitative  data  to  improve  the  growth  equation. 

The  second  specific  objective  involving  the  nucleation  threshold 
conditions  is  equally  important  in  view  of  the  catastrophic  nature  of 
the  adiabatic  shear  banding  process.  That  is,  in  armor  penetration 
scenarios  the  onset  of  shear  banding  in  a  significant  portion  of  the 
armor  is  often  tantamount  to  failure  of  the  armor  by  plugging.  Thus, 
an  armor  designer  must  attempt  to  prevent  a  specific  threat  from 
exceeding  the  threshold  conditions.  These  threshold  conditions  are  then 
the  governing  conditions. 

The  third  specific  objective,  that  of  establishing  optimal  experi¬ 
mental  and  analytical  procedures  for  measuring  the  key  material  properties 
governing  shear  band  behavior,  is  important  for  evaluating  and  character¬ 
izing  candidate  armor  and  penetrator  materials. 

1  2 

The  SNAG  shear  band  model  developed  previously  ’  was  based  on 
observations  of  shear  banding  stopped  at  early  stages  in  contained 
fragmenting  cylinder  experiments  similar  to  those  described  in  detail 
in  the  next  section. 

The  experiments  performed  on  HF-1  high  fragmentation  steel  indicated 
that  shear  bands  nucleate  and  grow  as  a  function  of  the  plastic  strain 
and  strain  rate.  Bands  usually  started  at  the  inner  diameter  of  the 
round,  proceeded  outwards  at  45  degrees  to  the  radial  direction,  and 
stayed  approximately  semicircular.  The  shear  displacement  in  a  band 
was  proportional  to  the  band  length.  Small  ripples  or  striations 
nucleated  in  HF-1  at  15  percent  plastic  strain  and  began  to  grow  at 
about  30  percent  plastic  strain. 
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The  computational  model  constructed  to  represent  the  experimental 
phenomena  contains  nucleation  and  growth  of  shear  bands,  plus  a  damage- 
caused  stress  reduction  process.  Nucleation  occurs  according  to  the  equation 


=  c  fU 

dt  n  n\E 


2  ,  PS 

dc , . 


where  N^q  is  the  number  of  shear  bands  per  unit  volume  on  the  <}>0  plane, 

t  is  time,  C  is  a  coefficient  with  the  units  of  time  squared,  F  is  a 
n  n  n 

factor  representing  the  fraction  of  solid  angle  assigned  to  the  00 

p 

direction,  E  is  a  critical  energy,  dE  /dt  is  the  rate  of  increase  of 
cr  PS 

plastic  strain  energy,  and  is  the  plastic  shear  strain  on  the  00 

plane.  The  shear  bands  are  not  of  uniform  size  but  are  nucleated  with 
an  exponential  size  distribution. 

Growth  occurs  by  the  viscous-type  growth  model : 

PS 

dF  ~  Cg  dt  R  (2) 


where  C  is  a  dimensionless  growth  coefficient,  and  R  is  the  nominal 
radius  of  a  band. 

3 

The  damage  on  any  plane  is  taken  as  proportional  to  ENR  on  that 
plane.  The  stress  that  can  be  sustained  on  a  plane  is  reduced  in  direct 
proportion  to  this  damage  factor. 

At  the  conclusion  of  the  simulation,  each  computational  cell  that 
has  undergone  shear  banding  contains  a  distribution  of  sizes  of  shear 
bands  on  one  or  more  planes.  These  distributions  are  assembled  to 
produce  a  fragment  size  distribution  for  the  cell  and  for  the  entire 
object . 

With  the  shear  band  model,  simulations  were  made  of  partially  and 
completely  fragmenting  cylinders  and  a  military  round.  The  correlation 
of  the  measured  fragment  size  distributions  with  the  computed  ones  was 
quite  good  for  test  conditions  that  included  two  sizes  of  cylinders, 
two  loading  rates  (confined  or  free),  and  two  materials  (HF-1  steel  and 
Armco  iron)  with  very  different  yield  strengths. 


At  the  conclusion  of  the  previous  shear  band  work,  it  was  recognized 
that  there  remained  several  gaps  in  our  knowledge  of  shear  banding.  The 
threshold  conditions  of  strain  rate,  strain,  and  stress  for  initiating 
shear  bands  were  not  understood.  Critical  temperatures  or  energies 
used  as  nucleation  parameters  were  simply  estimated.  A  growth  process 
was  assumed  that  would  be  consistent  with  the  final  size  distributions 
of  shear  bands,  but  there  was  no  experimental  basis  for  uniquely 
determining  the  growth  rate.  The  interaction  of  the  three  observed  shear 
band  modes  was  not  understood,  although  it  was  assumed  that  together 
the  three  modes  determine  the  fragment  size  and  shape.  Stress-strain 
relations  for  partially  shear-banded  material  were  derived  on  the  basis 
of  solid  mechanics  relations  plus  reasonable  physical  assumptions. 

However,  these  stress-strain  relations  should  receive  independent  experi¬ 
mental  verification  before  they  can  be  assumed  to  represent  the  actual 
damage  processes. 

To  attack  the  above  deficiencies  in  the  SNAG  model,  in  the  present 
program  we  have  concentrated  on  designing,  performing,  and  analyzing 
experiments  that  would  yield  optimal  information  on  shear  band  initiation 
and  growth.  As  discussed  later  in  Section  6,  contained  fragmenting  cylinder 
experiments  were  chosen  as  the  most  efficient  means  for  generating  the 
desired  data.  A  theoretical  model  of  shear  band  initiation  was  developed 
(Section  5);  this  theory  predicts  that  the  shape  of  the  work  hardening 
curve  for  the  material  is  the  governing  material  property  for  shear  band 
initiation.  Therefore,  4340  steels  of  differing  hardness  levels  and 
correspondingly  differing  work  hardening  curves  were  chosen  for  the 
experiments . 

By  varying  the  explosive  density  and  the  experimental  configurations, 
we  were  able  to  vary  the  stress  and  strain  histories  experienced  by  the 
material  elements  in  the  exploding  cylinders.  Computer  calculations  of 
these  load  histories  were  compared  with  observed  shear  band  damage, 
thus  allowing  initial  estimates  to  be  made  of  the  initiation  conditions 
and  the  shear  band  growth  velocities. 


In  the  remainder  of  this  report  we  describe  in  detail  the  progress 
we  have  made  toward  the  program  objectives.  Section  3  describes  the 
contained  fragmenting  cylinder  experiments,  and  presents  the  reduced 
data  from  these  experiments.  Section  4  discusses  the  computational 
simulations  of  the  experiments,  and  presents  our  first  estimates  of  the 
parameters  obtained  for  our  shear  banding  model,  and  Section  5  discusses 
our  proposed  shear  band  nucleation  criterion.  Section  6  discuss  optimal 
experimental  test  configurations  for  characterizing  the  shear  band 
behavior  of  materials,  and  in  Section  7  we  discuss  the  results  and 
present  our  recommendations  for  further  work. 
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3  CONTAINED  FRAGMENTING  CYLINDER  EXPERIMENTS 


Experimental  Technique 

To  obtain  quantitative  shear  banding  data  on  4340  steel,  we  used  the 
contained  fragmenting  cylinder  experimental  technique  devised  under  a 
previous  program^  for  study  of  shear  banding  failure  in  cylindrical  pipes. 
This  technique  enables  us  to  load  thick-walled  cylindrical  tubes  at  pres¬ 
sures  and  strain  rates  high  enough  to  induce  shear  banding,  but  to  contain 
the  tubes  in  such  a  manner  as  to  arrest  the  shear  banding  process  at  a 
level  before  full  fragmentation  occurs.  We  are  then  able  to  examine  the 
material  metallographically ,  observe  the  shear  bands,  and  relate  them  to 
the  loading  histories. 

The  experimental  geometry  is  shown  in  Figure  1.  The  4340  steel 

* 

specimen  tube  is  filled  with  the  high  explosive  PETN  ,  pressed  in  1.27  cm 

3 

axial  increments  to  a  density  between  0.7  and  1.47  gm/cm  .  A  gap  is  left 
between  the  PETN  and  the  upper  and  lower  axial  end  of  the  specimen  tube 
in  order  to  prevent  unwanted  edge  effects  from  shearing  off  the  end  of 
the  tube.  Detonation  occurs  through  a  PETN  cone  to  produce  a 
relatively  planar  detonation  wave  traveling  down  the  tube  in  the  axial 
direction.  The  high  internal  pressures  from  the  PETN  deform  the  tube 
outward,  inducing  a  state  of  strain  that  has  a  large  compressive  radial 
component,  a  large  tensile  hoop  component,  and  a  negligible  axial  component 
(ignoring  end  effects).  We  thus  have  a  state  of  high  shear  strain. 

The  specimen  tube  is  surrounded  by  a  series  of  concentric  contain¬ 
ment  pipes.  The  relatively  thin-walled,  soft  clear  acrylic  buffer  tube 
allow  the  specimen  cylinder  to  expand  radially  and  prevents  a  sharp  high- 
pressure  impact  with  the  main  containment  tube,  which  is  a  very  thick- 
walled  mild  steel  pipe.  The  outermost  lead  pipe  serves  as  a  momentum 


* 


Pentaerythritol  tetranitrate  (C^HgO-j  2^ 


). 
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FIGURE  1  CROSS-SECTION  OF  CONTAINED  FRAGMENTING  CYLINDER  EXPERIMENTS 
X  varies  from  0.64  to  1.9  cm,  p  varies  from  0.7  to  1.47  gm/cm3. 


trap;  it  separates  from  the  mild  steel  pipe  and  moves  out  in  the  radial 
direction,  taking  most  of  the  radial  momentum  with  it  and  reducing  the 
magnitude  of  the  inward-going  tensile  wave.  Because  of  the  thickness  of 
the  containment  tubes,  the  specimen  tube  experiences  no  radial  tensile 
stress  for  at  least  40  psec,  by  which  time  the  internal  pressures  produced 
by  the  PETN  are  greatly  reduced,  due  to  escape  of  explosive  gases  from 
the  open  ends  of  the  tube.  The  lack  of  radial  tensile  stresses  inhibits 
brittle  fracture,  leaving  shear  banding  as  the  predominant  failure  mode 
in  the  specimen  tube. 

The  strain  and  strain  rate  to  which  the  specimen  tube  is  subjected 
are  varied  by  changing  the  explosive  mixture  and  the  thickness  of  the 
acrylic  buffer  tube.  A  higher  pressure  detonation  will  induce  higher 
strain  rates  in  the  specimen  tubes,  and  a  thicker  buffer  will  allow  the 
specimen  tube  to  attain  larger  strains.  However,  the  controlling 
parameters  are  not  orthogonal  with  respect  to  the  material  response 
parameters.  For  example,  a  higher  detonation  pressure  will  produce 
larger  strains  even  if  the  buffer  tube  thickness  remains  the  same. 
Furthermore,  there  is  a  complicated  pattern  of  stress  waves  and  reflec¬ 
tions  while  the  various  containment  tubes  are  "ringing  up"  to  their 
peak  pressures  and  then  unloading.  Hence,  the  loading  history  of  the 
specimen  tube  cannot  be  predicted  easily  without  the  aid  of  computer 
simulations.  In  section  4  we  present  the  results  of  the  computer 
simulations  and  show,  for  example,  how  the  computed  strain  rates  vary 
with  detonation  pressure. 

A  total  of  seventeen  contained  fragmenting  cylinder  experiments 

were  performed  in  this  program,  using  commercially  obtained  tubular  4340 

steel  as  the  specimen  material.  Eight  shots  were  fired  using  the  4340 

steel  treated  to  R  40,  seven  used  steel  treated  to  R  52,  and  two  used 
c  c 

fully  annealed  steel  (~Rc21).  The  experimental  parameters  for  each  shot 
are  given  in  Table  1  and  are  also  depicted  graphically  in  Figure  2, 
which  shows  the  location  of  each  experiment  in  PETN  density — acrylic 
buffer  thickness  space.  On  this  graph,  a  positive  change  in  the  hori¬ 
zontal  position  would  represent  a  greater  shear  strain  loading  rate  in 
the  specimen,  and  a  positive  vertical  change  would  represent  a  greater 
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loading  duration.  Either  change  or  a  combination  of  both  should  lead  to 
a  greater  level  of  shear  band  damage. 

The  loading  parameters  for  the  Rc40  and  R^52  experiments  were 
selected  to  cover  a  broad  range  of  damage  from  negligible  shear  banding 
(homogeneous  deformation  only)  to  substantial  fragmentation.  Qualita¬ 
tive  damage  results  are  given  in  Table  1.  The  experimental  series  were 
also  designed  so  that,  for  several  cases,  experiments  having  identical 
loading  parameters  would  be  performed  using  specimen  materials  of  the 
different  hardnesses,  so  as  to  compare  the  resistance  to  shear  banding 
as  a  function  of  material  hardness  or  yield  strength. 

Qualitative  Results 

For  the  shots  resulting  in  less  damage,  the  specimen  tube  was 
recovered  in  one  or  two  large  pieces,  which  we  then  sliced  parallel  tc 
the  axis  to  make  the  inner  surface  more  visually  accessible  (Figures  3a 
and  3b).  For  the  shots  resulting  in  more  damage,  the  tube  was  fragmented 
into  as  many  as  30  pieces.  The  recovered  fragments  were  arranged  into 
a  jig-saw-puzzle-type  pattern  according  to  their  relative  position  in 
the  original  specimen  tube  (e.g.,  as  in  Figure  3c),  so  that  we  could 
determine  the  damage  as  a  function  of  axial  position  in  the  tube.  Com¬ 
parison  of  the  fragments  recovered  from  shots  with  identical  loading 
parameters,  but  with  specimen  materials  of  different  hardness  (as  in 
Figure  3),  showed  that  the  harder  material  will  shear  band  more  readily. 
The  fully  annealed  (R^21)  4340  steel  exhibited  no  shear  banding  (except 
for  three  large  shear  cracks  that  originated  at  or  near  the  upper  end 
of  the  tube  as  a  result  of  edge  effects  altering  the  basic  strain  pattern) 
despite  having  undergone  a  very  large  amount  of  homogeneous  plastic 
deformation.  In  contrast,  the  substantially  hardened  R^40  steel 
exhibited  numerous  small  shear  bands  widely  scattered  over  most  of  the 
inner  surface,  although  it  experienced  only  about  half  as  much  homo¬ 
geneous  plastic  deformation,  and  the  extremely  hardened  R  52  steel 

c 

exhibited  extensive  shear  banding  leading  to  fragmentation  after  even 
less  homogeneous  plastic  deformation.  In  Section  5,  we  discuss  reasons 
for  this  effect. 
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The  three  orientation  modes  for  shear  bands  in  an  expanding  cylinder 
are  shown  in  Figure  4.  As  discussed  in  detail  by  Erlich  et  al.,*"  the 
modes  represent  shear  slip  along  the  plane  of  maximum  shear  strain.  In 
this  present  study  of  4340  steel,  as  well  as  in  the  previous  study  of 
HF-1  steel,  the  mode  1  shear  bands  predominate.  They  are  the  only  ones 
evident  in  the  experiments  producing  lower  damage  levels.  At  the  higher 
levels  of  damage,  although  mode  2  bands  appear  at  the  inner  surface  in 
significant  numbers  and  mode  3  bands  are  occasionally  seen  in  the 
interior,  it  is  still  the  mode  1  bands  that  have  become  the  longest  and 
have  taken  up  the  most  strain.  The  mode  1  bands  extend  all  the  way 
through  a  specimen  or  coalesce  with  other  mode  1  bands  to  cause  frag¬ 
mentation  into  axial  strips,  although  the  mode  2  bands  may  determine 
the  length  or  aspect  ratio  of  these  strips. 

In  the  previous  shear  banding  study  of  HF-1  steel,’*'  it  was  determined 
that  more  than  two-thirds  of  the  mode  1  shear  bands  created  in  a  contained 
fragmenting  cylinder  experiment  intersected  the  inner  surface  of  the 
specimen  tube.  It  was  also  determined  that  the  length  ( £)  of  the  shear 
band  in  the  direction  perpendicular  to  its  slip  plane  (see  Figure  4) 
was  approximately  twice  its  depth  (d)  along  the  slip  plane  and  that  the 
depth  was  roughly  15  times  the  maximum  shear  displacement  (B)  along  the 
slip  plane. 

To  determine  if  these  findings  also  hold  true  for  4340  steel,  we 
selected  a  fragment  from  shot  2,  shown  in  Figure  5,  and  made  a  series  of 
seven  cuts  perpendicular  to  the  tube  axis  with  a  thin  diamond  blade.  We 
polished  and  etched  each  cut  surface,  identified  the  mode  1  shear  bands 
common  to  adjacent  cuts  and  measured  the  shear  band  depth  and  displace¬ 
ment  on  each  cut.  The  results  for  one  shear  band,  labeled  no.  6,  are 
given  in  Figure  5c,  which  shows  d  and  B  plotted  against  the  relative 
distance  along  The  band  appears  to  have  a  half-length  of  almost 
6-1/2  mm,  a  maximum  depth  of  about  7-1/2  mm,  and  a  maximum  shear  dis¬ 
placement  of  ~  2/3  mm.  Its  shape  is  very  nearly  that  of  a  semicircle; 
thus,  if  a  band  is  assumed  to  nucleate  at  one  point,  it  grows  in  depth 
about  as  much  as  it  grows  in  length  in  each  direction.  The  depth  appears 
to  be  on  the  order  of  11  or  12  times  the  displacement.  In  addition,  all 
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FIGURE  4  GEOMETRY  AND  NOMENCLATURE  FOR  SHEAR  BANDS  IN  FRAGMENTING 
CYLINDER  EXPERIMENTS 


SHEAR  BAND  DIMENSION 


FIGURE  5  FRAGMENT  (a)  AND  (b>  USED  IN  SHEAR  BAND  PROFILE  STUDY,  AND 
MEASURED  PROFILE  (c)  OF  SHEAR  BAND  NO.  6 


the  mode  1  bands  seen  on  the  cuts  in  this  fragment  or  any  other  4340 
steel  fragment  intersected  the  inner  surface  of  the  specimen  tube. 

It  was  noticed  that  for  shots  resulting  in  small  amounts  of  shear 
banding  damage,  the  specimen  tube  would  sometimes  separate  along  a 
diagonal  line  that  would  around  the  tube  at  approximately  15°  to  20° 
from  the  axial  direction  (as  in  Figure  3b).  For  shots  with  more  sub¬ 
stantial  shear  banding  damage,  often  the  shear  bands  would  link  up  along 
this  same  diagonal  direction,  although  each  individual  band  was  formed 
parallel  to  the  axis.  It  is  believed  that  this  effect  is  caused  by 
texturing  of  the  4340  steel  tube  in  the  process  of  extrusion. 

Quantitative  Results 

As  a  result  of  the  findings  discussed  above,  we  wore  able  to 
simplify  our  method  for  determining  quantitative  shear  band  damage. 

Rather  than  making  a  large  number  of  cuts  in  each  fragment  and  polishing, 
etching,  and  photographing  the  surface  of  each  cut,  and  then  measuring 
the  depth  of  each  shear  band,  we  decided  to  photograph  in  detail  the 
inside  surface  of  each  fragment  and  measure  the  length  and  axial  position 
of  each  shear  band.  This  procedure  is  much  less  time  consuming. 

As  an  example,  Figure  6  shows  photomicrographs  of  part  of  one 
fragment  recovered  from  Shot  5  that  were  used  in  the  quantitative  shear 
band  determination.  Two  photomicrographs  with  side  lighting  from 
opposite  directions  (Figure  6b  and  6c)  are  used  to  highlight  the  mode  1 
shear  bands  intersecting  the  surface  at  the  two  orthogonal  45  degree 
angles.  The  fragment  is  divided  into  zones  of  axial  position,  each  zone 
being  1.27  cm  wide,  and  shear  bands  within  each  zone  are  counted  as  a 
function  of  length.  Values  for  each  zone  are  added  together  for  all  of 
the  fragments  recovered  for  a  particular  lot,  and  the  total  is  divided 
by  the  surface  area  of  the  zone  to  yield  shear  band  surface  densities. 

Quantitative  damage  assessments  were  made  only  for  the  shots  using 
4340  steel  specimen  tubes  hardened  to  Rc40  (except  for  Shots  1  and  6, 
which  exhibited  negligible  shear  banding).  The  results  are  shown  in 
Figures  7  through  12,  plotted  in  terms  of  the  areal  density  of  shear 
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FIGURE  6  FRAGMENT  (a)  RECOVERED  FROM  SHOT  NO  5.  AND  PHOTOMICROGRAPHS  (t>  and  r) 
OF  TWO  ZONES  FOR  USE  IN  QUANTITATIVE  SHEAR  BAND  STUDY 
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FIGURE  8  SHEAR  BAND  DISTRIBUTION  FOR  SHOT  3 
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SURFACE  DENSITY  OF  BANDS  WITH  LENGTH  >  L  —  cm 
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FIGURE  9  SHEAR  BAND  DISTRIBUTION  FOR  SHOT  4 
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SURFACE  DENSITY  OF  BANDS  WITH  LENGTH  >  L  —  cm': 
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FIGURE  12  SHEAR  BAND  DISTRIBUTION  FOR  SHOT  8 
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bands  having  a  length  greater  than  £  versus  the  shear  band  length 
for  different  zones  of  axial  position. 

In  examining  the  graphs,  one  should  note  several  features.  First, 
the  shear  band  densities  are  higher  for  zones  near  the  axial  mid-point 
of  the  specimen  tube  than  for  zones  near  either  end.  This  is  because  the 
pressures  from  the  explosives  remained  higher  for  a  longer  duration. 
Second,  the  very  small  shear  bands  having  a  length  less  than  about  0.2  cm 
are  difficult  to  see,  since  the  width  of  their  exposed  slip  plane  (B, 
from  Figure  5)  may  be  less  than  0.1  mm;  hence,  the  number  of  these  bands 
counted  may  be  below  the  actual  value.  Therefore,  the  shear  band  density 
curves  for  values  of  £  less  than  0.2  cm  (see  Figure  7,  for  example)  may 
actually  be  as  steep  as  they  are  for  larger  values  of  £ 

In  summary,  we  have  performed  a  series  of  contained  fragmenting 
round  experiments  on  specimens  of  4340  steel  treated  to  three  different 
hardnesses  and  have  observed  a  range  of  damage  from  homogeneous  plastic 
deformation  only  to  extensive  shear  banding  and  fragmentation.  We 
devised  a  simplified  method  of  quantifying  the  shear  banding  damage 
and  obtained  shear  band  density  distributions  as  a  function  of  axial 
position  for  each  of  the  Rc40  steel  experiments.  Quantification  of  the 
shear  band  damage  for  the  1*^52  steel  experiments  is  recommended  for 
future  work. 
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4  COMPUTATIONAL  ANALYSIS 


The  TROTT  computer  code,  which  was  developed  recently  at  SRI  and 
3 

documented  by  Seaman,  was  used  to  simulate  all  the  R  40  contained 

c 

fragmenting  round  experiments  for  which  we  had  shear  band  density  dis¬ 
tributions.  The  TROTT  code  is  a  two-dimensional  Lagrangian  wave  propaga¬ 
tion  code  incorporating  most  of  the  complex  equations-of-state  and 

4 

fracture  models  available  in  the  one-dimensional  PUFF  code  .  The 
original  calculations  with  the  TROTT  code  did  not  utilize  a  slide  line 
at  the  HE-specimen  tube  interface.  Therefore,  in  the  simulations  the 
explosive  gases  are  not  allowed  to  escape  from  the  top  and  bottom  of 
the  specimen  tube  as  fast  as  they  would  in  the  actual  experiment. 

However,  this  was  not  expected  to  greatly  affect  the  pressure  histories 
in  the  mid-axial  regions  until  reasonably  late  in  the  simulation  (=30 
psec),  so  that  the  response  of  the  specimen  tube  in  the  mid-axial  zones 
(zones  D  &  E  in  Figures  7  through  12)  may  be  approximated  fairly  well 
by  the  TROTT  code  without  the  slide  line. 


The  equation  of  state  for  reduced  density  PETN  was  obtained  from 

the  experimental  data  of  Hornig,  Lee,  Finger,  and  Kurrle.^  They  measured 

Chapman-Jouguet  pressures  (peak  detonation  pressures)  and  detonation 

3  3 

velocities  in  PETN  for  densities  from  0.25  g/cm  to  1.764  g/cm 
(approximately  the  crystal  density).  For  our  calculations  the  data 
were  fitted  to  a  polytropic  gas  form  to  obtain  the  polytropic  gas 
constant  y  and  the  chemical  energy  Q. 


To  fit  the  data  to  the  polytropic  form,  the  detonation  velocity  D 
and  pressure  PCJ  were  plotted  versus  density  as  shown  in  Figure  13.  A 
smooth  curve  was  drawn  through  the  points.  The  values  of  y  and  Q  were 
obtained  from  the  polytropic  gas  relations 


Y  = 


P  D 
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-  1 


and  Q  = 
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FIGURE  13  DETONATION  VELOCITY  (a)  AND  CHAPMAN-JOUGUET  PRESSURE 

<b|  FOR  PETN  AS  A  FUNCTION  OF  DENSITY, 

ACCORDING  TO  HORNIG  et  al.5 


I 


where  D  =  detonation  velocity,  cm/sec 

2 

P  =  Chapman-Jouguet  pressure,  dyn/cm 

Q  =  chemical  energy,  erg/g 

3 

=  initial  density,  g/cm 

The  resulting  constants  for  the  PETN  are  shown  in  Figure  14.  These 
constants  were  used  in  our  simulations  of  fragmenting  round  experiments. 

The  simulations  were  made  using  the  shear  banding  subroutine 
SHEAR2^  so  as  to  calculate  the  plastic  shear  strain  in  directions 
appropriate  to  each  of  the  three  orientation  modes  for  shear  banding  in 
a  cylindrically  expanding  pipe.  The  parameters  were  set  in  these  calcu¬ 
lations  so  as  not  to  allow  any  shear  bands  to  nucleate  or  grow.  Thus, 
the  simulations  were  for  the  case  of  negligible  shear  banding.  They 
would  be  expected  to  be  accurate  for  experiments  that  produced  negligible 
or  incipient  shear  banding  or  for  experiments  that  produced  higher  damage 
levels  before  the  onset  of  significant  band  growth,  at  which  time  sign¬ 
ificant  strains  would  be  taken  up  by  the  shear  bands  themselves.  In 
any  case,  the  simulations  would  give  a  first  approximation  to  the  strain 
histories . 

Since  the  shear  band  density  distributions  that  we  measured 
quantitatively  from  the  Rc40  experiments  involved  mode  1  shear  bands 
only,  we  calculated  mode  1  strain  histories  (Ref.  1,  p.  35).  The  strain 
rate  history  was  selected  because  the  strain  rate  is  the  parameter 
analogous  to  the  stress  in  the  case  of  tensile  failure.  Just  as  the 

6 

growth  rate  of  tensile  voids  and  cracks  depends  on  the  tensile  stress, 
we  expect  that  the  shear  band  growth  rate  depends  on  the  shear  strain 
rate . 

The  plastic  shear  strain  rate  histories  and  the  shear  band  density 
distributions  for  the  mid-axial  location  (an  average  of  zones  D  &  E  in 
Figures  7  through  12)  are  grouped  according  to  the  explosive  density 
in  Figures  15  through  17.  Figure  15  includes  Shots  3  and  5,  both  of 

3 

which  had  an  H.E.  density  of  1.35  g/cm  ;  Figure  16  includes  Shots  2,  4, 
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POLYTROPIC  GAS  CONSTANT, 


FIGURE  14  POLYTROPIC  GAS  CONSTANT  y  AND  CHEMICAL  ENERGY  Q  REDUCED 

FROM  DATA  OF  HORNIG  et  al.5  ON  PETN  AT  SEVERAL  DENSITIES 
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FIGURE  15  COMPUTED  MID-AXIAL  MODE  1  PLASTIC  SHEAR  STRAIN  RATE  HISTORIES 
la)  FOR  SHOTS  3  AND  5,  COMPARED  WITH  EXPERIMENTALLY  DETERMINED 
SHEAR  BAND  DENSITY  DISTRIBUTIONS  (b). 


FIGURE  16  COMPUTED  MID-AXIAL  MODE  1  PLASTIC  SHEAR  STRAIN  RATE  HISTORIES 

(a)  FOR  SHOTS  2,  4,  AND  7.  COMPARED  WITH  EXPERIMENTALLY  DETERMINED 
SHEAR  BAND  DENSITY  DISTRIBUTIONS  (b). 


SCJNV8  30  A1ISN3Q  33V3HnS 


MS/frOl  —  31VU  NIVUIS 


O 

00  Z 

r** 

1-  2 

o 

E 

§ 

o  t- 
X  3 

o 

I 

? 

< 

(0  m 

X  X 
O  H 

00 

c* 

5 

LL  tO 

o 

*5  ° 

X* 

~  > 

CO 

t— 

C3 

>  H 
X  Jo 
O  Z 

o 

Z 

t-  UJ 

LLI 

■Q 

to  Q 

— 1 

1  Q 

Q 

UJ  Z 

o 

Z 

< 

I-  < 
<  m 
nr 

CQ 

X 

2  < 

CM 

DC 

~  UJ 

o 

< 

LLI 

%  X 
X  to 

I 

to  Q 

C/) 

„  ul 

o 

<  ? 

uj  S 

I  X 
to  Lu 

2  t 

<  > 
_i  _i 

0 - 1 

_  < 

kO 

H 

CN 

UJ  Z 
□  UJ 

O  5 

O 

S  x 

CM 

a 

m 

_l  UJ 

<  °- 
—  X 

IT) 

X  uj 

1 

(a) 

*?  X 
D  1- 

O 

LLI 

2 

i  5 

Q  o 

1- 

UJ  UJ 
1-  X 

IT) 

3  < 
a.  cl 
5  5 

O  O 

O 

V  o 

r- 

• — 

LU 

X 

3 

O 

39 


_ 
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and  7,  all  of  which  had  an  H.E.  density  of  1.2  g/cm  ;  and  Figure  17 

3 

depicts  Shot  8,  which  had  an  H.E.  density  of  1.47  g/cm  . 

The  peak  strain  rates  are  almost  exactly  the  same  for  shots  having 

the  same  H.E.  density  and  hence  the  same  detonation  pressures.  These 

peak  strain  rates  are  shown  in  Figure  18  as  a  function  of  detonation 

pressure.  The  rates  range  from  0.7  x  10^/sec  for  a  detonation  pressure 

4 

of  33  kbar  to  2.2  x  10  /sec  for  a  detonation  pressure  of  204  kbar . 

For  shots  of  a  given  H.E.  density,  increasing  the  acrylic  buffer  thickness 
increases  the  shear  strain  rate  duration  but  not  the  peak  amplitude. 

The  change  in  the  computed  pulse  duration  appears  to  vary  linearly  with 
the  increase  in  buffer  thickness  and,  in  fact,  has  a  value  equal  to  the 
time  it  takes  a  stress  pulse  to  traverse  back  and  forth  across  the 
change  in  thickness. 

From  the  shear  band  density  distributions  and  the  strain  rate 
histories  shown  in  Figure  15  and  16,  it  appears  that: 

(1)  Higher  strain  rates  result  in  greater  shear  band 
nucleation  rates. 

(2)  A  slight  increase  in  the  duration  of  the  peak  strain 
rates  causes  a  marked  increase  in  the  shear  band  lengths. 

(3)  The  shear  band  growth  rates  are  larger  for  higher  peak 
strain  rates. 

We  noticed  that  the  shear  band  density  distribution  curves  are 
approximately  straight  lines  in  the  semi-log  graph.  So  in  our  first 
attempt  to  quantify  these  qualitative  results,  and  to  obtain  parameters 
that  would  be  used  in  the  initial  computer  simulation  in  which  shear 
banding  was  allowed,  we  fitted  the  data  with  straight  lines  such  that 
lines  for  experiments  with  the  same  peak  strain  rates  have  the  same 
y-axis  intercept.  The  results  are  shown  graphically  in  Figure  19.  A 
consequence  of  this  procedure  is  that  the  shear  bands  follow  a  viscous 
growth  law  of  the  form 


R/R  =  constant. 
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SHEAR  BAND  LENGTH,  C  —  cm 

MA-6408-24 

FIGURE  19  SHEAR  BAND  DENSITY  DISTRIBUTIONS  FORCED  TO 
OBEY  VISCOUS  GROWTH  LAW 
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where  R  is  the  shear  band  radius  (equal  to  i$L)  .  This  appears  at  first 
glance  to  be  in  agreement  with  the  growth  law  assumed  in  References 
1  and  2  and  given  by  Eq.  (2)  in  Section  2.  However,  the  initial  evidence 
suggests  that  R/R  is  not  a  linear  function  of  strain  rate,  since  for  the 

^  pc  4 

cases  shown  in  Figure  19,  C  S'  21  for  e,n  =  2.05  x  10  /sec,  whereas  C 

~  PS  48  (J>e  g 

=  8  for  =  1.7  x  10  /sec.  Subsequent  iterative  calculations  may  show 

that  a  constant  value  of  C  can  describe  the  data,  but  at  present  it 

•  g 

appears  R/R  will  be  a  more  complicated  function  of  strain  rate  than  that 
given  by  Eq.  (2). 

The  magnitude  of  the  shear  band  growth  velocity  approaches  1.7  mm/ys 
for  the  highest  strain  rates.  This  speed  is  about  half  the  shear  wave 
velocity  for  4340  steel  («3.2  mm/ys),  which  is  the  theoretical  limit  for 
such  instabilities. 

As  expected,  the  nucleation  rate  appears  to  be  a  strong  function 

of  strain  rate.  The  total  surface  density  of  nucleated  shear  bands 

-2  -2 

(see  Figure  19)  increases  from  16  cm  to  125  cm  as  the  nominal  peak 

4  4 

strain  rate  increases  from  1.7  x  10  /sec  to  2.05  x  10  /sec. 

At  the  very  end  of  the  program,  after  a  slide  line  capability  had 
been  inserted  into  the  TROTT  code,  we  re-ran  the  no-damage  calculation 
of  Shot  2,  using  a  slide  line  to  simulate  the  motion  of  the  explosive 
gases  along  the  interface  with  the  specimen  tube.  To  illustrate  this 
calculation,  we  plotted  various  motions,  stresses  and  strains  experienced 
by  the  specimen  tube  at  its  mid-axial  location.  The  cell  layout  is  shown 
in  Figure  20,  along  with  the  location  of  the  cells  for  which  we  have 
plotted  the  histories.  The  HE  begins  detonation  simultaneously  across 

the  top  and  the  detonation  travels  downward.  Figure  21  shows  the 

* 

pressure  history  in  the  HE  cell  adjacent  to  the  specimen  tube  near  the 
mid-axial  location.  The  radial  position  and  velocity  histories  of  the 
4  mid-axial  cells  (J  =  5  through  8)  are  shown  in  Figure  22.  The  principal 

•k 

stress  histories  within  the  specimen  tube  (cells  J  =  6  through  8)  are 


Note  that  positive  values  for  pressures  and  stresses  connote  compression; 
negative  values  connote  tension. 
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LEAD  MILD  STEEL  LUCITE  \ 

MOMENTUM  CONTAINMENT  BUFFER  SPECIMEN 


TRAP  TUBE  TUBE  TUBE 


Location  of  Histories  Plotted  in  Subsequent  Figures  for: 

•  Positions  and  Velocities,  a  Stresses  and  Strains,  *  Pressures 

M A -6408-30 


FIGURE  20  CELL  LAYOUT  FOR  COMPUTER  SIMULATION  OF  CONTAINED  FRAGMENTING 
CYLINDER  EXPERIMENTS 
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PRESSURE  —  kbar 


FIGURE  21  CALCULATED  PRESSURE  HISTORY  IN  THE  H.E.  ADJACENT  TO  THE 
SPECIMEN  TUBE  NEAR  THE  MID-AXIAL  LOCATION 
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VELOCITY  —  mm/MSec  DISTANCE  FROM  AXIS  —  cm 


FIGURE  22  CALCULATED  RADIAL  POSITION  (a)  AND  VELOCITY  (b)  HISTORIES 
FOR  SPECIMEN  TUBE  AT  THE  MID-AXIAL  LOCATION 
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are  shown  in  Figure  23.  The  equivalent  plastic  strain*  histories  and 
the  plastic  shear  strain  histories  in  the  mode  1,  2,  and  3  directions 
(as  defined  in  Figure  4)  are  shown  in  Figures  24  and  25,  respectively. 

A  comparison  of  this  slide-lined  calculation  with  the  one  described 
earlier  for  Shot  2,  which  lacked  a  slide  line,  showed  that  the  slide- 
lined  calculation  exhibited  a  decrease  of  10%  in  a  peak  Mode  1  shear 
strain  rate,  a  decrease  of  18%  in  the  final  Mode  1  shear  strain,  and 
an  increase  of  4%  in  the  final  specimen  tube  thickness.  Furthermore, 
the  final  tube  thickness  attained  at  the  mid-axial  location  in  the 
slide-lined  calculation  was  within  1%  of  that  measured  experimentally. 
Thus,  the  original  calculations  without  the  slide  line  appear  to  be 
approximately  valid,  as  hoped.  However,  the  calculation  with  a  slide 
line  represents  a  meaningful  modification  to  the  contained  fragmenting 
cylinder  simulation  and  should  be  used  as  a  starting  point  for  future 
iterative  calculations  with  damage  that  will  be  necessary  to  better 
define  the  shear  banding  nucleation  and  growth  parameters  for  4340 
steel . 


* 

The  equivalent  plastic  strain  is  defined  as  3/2  times  the  square 
root  of  the  sums  of  the  squares  of  the  plastic  strains  in  the  three 
principal  directions. 
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MA-64 08-33 


FIGURE  23 


CALCULATED  PRINCIPAL  STRESS  HISTORIES  IN  THE 
SPECIMEN  TUBE  AT  THE  MID-AXIAL  LOCATION 


TIME  —  ix  sec 


M  A -6408-34 


FIGURE  24  CALCULATED  EQUIVALENT  PLASTIC  STRAIN  HISTORIES 
IN  THE  SPECIMEN  TUBE  AT  THE  MID-AXIAL  LOCATION 
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MA-64 08-35 


FIGURE  25  CALCULATED  PLAST'C  SHEAR  STRAIN  HISTORIES  IN 
THE  SPECIMEN  TUBE  AT  THE  MID-AXIAL  LOCATION 


5  SHEAR  BAND  NUCLEATION  THRESHOLD  CRITERION 


Introduction 

1  2 

In  our  original  development  of  the  SNAG  model  ’  we  did  not  attempt 
to  determine  what  causes  adiabatic  shear  bands  to  form  initially.  In¬ 
stead,  we  constructed  a  model  for  the  nucleation,  growth,  and  coalescence 
of  the  bands  once  they  appeared.  The  threshold  criterion  governing  the 
initial  appearance,  or  nucleation,  of  the  bands  was  estimated  empirically, 
and  it  stated  that  shear  bands  nucleate  when  both  a  critical  strain  and 
a  critical  strain  rate  are  exceeded.  No  derivation  was  presented. 

However,  such  a  derivation  of  a  nucleation  criterion  that  depends  on 
other,  commonly  measured  material  properties  would  be  very  desirable 
because  the  propensity  of  various  materials  to  form  shear  bands  could 
then  be  assessed  without  additional  extensive  testing. 

This  section  describes  our  recent  efforts  to  derive  and  test  such 
a  criterion.  We  first  give  some  background  information,  then  state  our 
postulated  criterion,  and  next  derive  the  mathematical  expressions  for 
the  nucleation  criterion.  Finally,  we  apply  the  criterion  to  several 
steels,  and  test  it  against  experimental  data. 


Background 


Classical  instability  analysis  is  based  on  assuming  the  existence 

7,8 

of  perturbations,  and  then  deriving  the  conditions  for  their  growth. 

The  usual  result  for  plasticity  problems  is  that  shear  instabilities  will 

occur  when  increasing  strain  results  in  decreasing  load,  that  is,  where 

the  work  hardening  coefficient  becomes  zero  or  negative.  However, 

instabilities  with  positive  work  hardening  coefficients  have  been  shown 

9  10 

to  be  possible  by  Rudnicki  and  Rice  and  by  Asaro  and  Rice.  A  recent 

11 

review  of  these  results  has  been  given  by  Rice. 
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The  above  instabilities  are  basically  static,  whereas  experimental 

evidence  suggests  that  adiabatic  shear  band  initiation  requires  high 

strain  rates.  However,  plastic  instabilities  can  also  be  viewed 
8  11 

dynamically.  ’  For  strains  equal  to  or  larger  than  the  critical  value, 

the  plastic  shear  wave  velocity  becomes  zero  or  imaginary.  Such  strains 

cannot  be  propagated  away  from  the  region  of  load  application,  and  the 

deformation  becomes  localized  as  a  growing  perturbation.  The  classical 

example  of  such  an  instability  is  that  first  investigated  independently 

during  the  Second  World  War  by  von  Karman,  Taylor,  White,  Le  Van  Griffis, 

12 

Rakhmatvolin,  and  Shapiro.  A  review  of  their  work  is  given  by  Davies. 

The  above  workers  found  that  for  wires  subjected  to  tensile  impact, 
strains  larger  than  that  corresponding  to  the  ultimate  strength  cannot 
be  propagated,  and  the  strain  is  localized  at  the  impact  end  of  the  wire, 
where  necking  occurs. 

It  is  clear  that  this  dynamic  instability  is  closely  related  to 
the  static  necking  instability.  The  usual  language  used  to  described 
the  necking  instability  is  that  for  strains  larger  than  those  at  ultimate 
tensile  strength  (UTS),  the  decrease  in  load  due  to  the  reduction  in 
area  exceeds  the  increase  in  load  due  to  work  hardening,  thus  causing 
decreasing  load  with  increasing  displacement,  and  fulfilling  the  con¬ 
dition  for  a  plastic  instability.  However,  because  the  Lagrangian 
plastic  bar  velocity,  /  1 _  9_o  \  i  (where  o  and  e  are  engineering  stress 

\P0  9e/ 

and  strain),  is  the  governing  wave  velocity  in  this  problem,  the  strain 
at  UTS  is  also  the  strain  for  which  the  wave  velocity  becomes  imaginary. 

For  static  or  dynamic  compression  of  a  rod,  however,  there  is  no 
corresponding  "barreling"  instability  because  the  force-displacement 
curve  always  has  a  positive  slope.  Nonetheless,  extremely  localized 
"barreling"  of  impact  rods  is  often  observed.  This  occurs  because  the 
bar  velocity  can  become  small  compared  to  the  impact  velocity,  even 
through  it  is  non-zero.  Thus,  the  larger  strains  do  not  propagate  very 
far  down  the  rod  during  the  duration  of  the  load.  This  phenomenon  is 
therefore  a  strain  localization  without  being  an  instability.  Neverthe¬ 
less,  the  dynamic  localization  process  is  the  same  in  both  cases: 
"trapping"  of  the  plastic  wave. 
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Postulated  Nucleation  Criterion 

We  postulate  that  adiabatic  shear  banding  is  a  dynamic  phenomenon 
whereby  plastic  strain  localization  at  a  perturbation  first  occurs  by 
"trapping  the  plastic  wave,"  and  is  then  followed  by  thermal  softening 
as  the  imposed  global  deformation  is  accomodated  solely  within  the 
perturbation,  thus  concentrating  the  plastic  work  and  the  associated 
heating.  Hence,  our  postulated  initiation  mechanism  is  not  restricted 
to  cases  of  zero  or  negative  work  hardening  (or  to  the  static  cases  with 
positive  work  hardening  discussed  by  Rice  and  coworkers)  but  is  a  purely 
dynamic  phenomenon  .  It  can,  however,  be  considered  as  an  extension  of 
the  usual  perturbation  approach  to  account  for  the  finite  rate  at  which 
plastic  flow  can  occur. 

We  thus  formulate  the  shear  band  initiation  criterion  as  follows: 
Plastic  shear  localization  leading  to  shear  banding  may  occur  at 
perturbations  when  the  rate  of  applied  shear  deformation  exceeds  the 
ability  of  the  plasti*  shear  wave  to  transmit  significant  deformation 
to  adjacent  material.  This  criterion  will  allow  dynamic  localization 
to  occur  when  there  is  positive  work  hardening  because  the  rate  of  applied 
deformation  is  non-zero.  For  quasistatic  loads,  the  rate  of  applied 
deformation  is  nearly  zero,  and  the  above  criterion  will  reduce  to  the 
usual  static  one  of  zero  or  negative  work  hardening. 

The  above  criterion  is  only  a  necessary  condition;  an  associated 
requirement  is  that  the  shear  stress-strain  curve  must  have  a  zero  or 
negative  second  derivative.  That  is,  the  plastic  wave  velocity  must 
remain  constant  or  decrease  with  increasing  strain.  Otherwise,  there 
could  be  a  "shocking-up"  effect  that  might  prevent  wave  trapping. 

Further,  the  wave  trapping  must  be  severe  enough  to  cause  significant 
thermal  softening.  Nevertheless,  we  shall  assume  for  the  present  that 
the  necessary  condition  is  the  primary  initiation  criterion. 
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Relation  of  Material  Properties  to  Flow  Field  Velocities 

The  above  criterion  can  be  quantified  by  considering  the  type  of 
deformation  shown  in  Figure  26.  We  assume  that  a  perturbation  in  an 
otherwise  smoothly  varying  strain  field  exists  in  the  x^,  x^  plane  and 
is  represented  by  the  area  A  in  that  plane.  Across  this  perturbation 
at  a  time  t  =  0,  there  is  a  discontinuous  jump  in  shear  strain.  For 
simplification,  we  assume  that  in  the  vicinity  of  the  perturbation,  all 
of  the  displacement  is  in  the  x^,  direction,  and  is  given  by 

u1  =  x1  -  h1  ,  (3) 


where 


X1  =  X1  (4) 

and  h^  =  Lagrangian  position  coordinate  . 

The  Lagrangian  shear  strain,  particle  velocity,  and  shear  strain  rate 
are  defined  by 


Y12  =  -  3u1/9h2 

v  =  Bu^/St  (3) 

*12  =  "  9vi/9h2 


Thus,  the  differential  displacement  may  be  written 

dul  =  8h^  dh2  +  9“T  dt  =  ‘  Y12dh2  +  vldt  *  (6) 
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The  Lagrangian  propagation  velocity  associated  with  a  constant  value  of 
u1  is  therefore 


Suj/at 

8u1/3h2 


(7) 


Similar  arguments  show  that  paths  of  constant  particle  velocity  v^  and 
shear  strain  y  have  the  velocities 


,  _  _  3v^/3t 

V1  3v1/3h2 


(8) 


and 


3Y12/9t 

Y12  9yl2/9h2 


'12 


9yi2/9h2 


(9) 


respectively. 

These  velocities  simply  describe  the  dynamic  deformation  field,  and 
do  not  contain  information  about  material  properties.  The  material 
properties  enter  through  a  constitutive  relation: 

T12  “  \2  <yij>  <10> 

where  is  the  shear  stress  in  the  h^  direction  across  the  perturba¬ 

tion  plane  represented  by  A  in  Figure  26.  Then  we  write  the  conservation 
of  momentum  equation  in  Lagrangian  coordinates: 

po(9vl/9t)h2  +  (^12/3h2)t  =  0. 
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In  the  vicinity  of  the  perturbation  A,  this  can  be  written  with  the  aid  of 
(N8)  as 


po(9vl/9t)h„ 


(3T12/^12)t(^12/8h2)t  =  °* 


(11) 


because  near  the  perturbation  all  strains  except  Y12  are  zero ’ 

Substituting  the  expressions  for  and  Y^2  from  (5)  into  (11), 
yields 


a2ul/3t2  2  Jj_  (12) 

32Ul/3h22  '  P  '  °» 


where  C 

P 


=  H/Po 


and  H  is  the  local  hardening  modulus.  An  identity  arising  from  the 
definitions  (5),  (8),  and  (9)  is: 


7  2 

3  u:/3t 

7  ? 

3  Ul/3h2^ 


(13) 


Comparison  of  (12)  and  (13)  yields  the  important  link  between  the 
deformation  field  and  the  work  hardening  curve  of  the  specific  material: 


(14) 


Examination  of  Flow  Field  for  Wave  Trapping 

To  define  the  conditions  for  wave  trapping,  we  assume  that  the 
perturbation  surface  represented  by  A  in  Figure  26  is  located  at 
h„  =  h  .  At  t  =  0  unstrained  material  on  one  side  of  the  perturbation 


on  the  other  side  of 


is  separated  from  material  at  a  finite  strain  y ^ 
the  perturbation.  We  further  assume  that  the  constitutive  relations  are 
the  same  for  the  material  on  each  side  of  the  perturbation. 

This  strain  discontinuity  introduced  at  t  =  0  will  cause  waves  to 
propagate  away  from  the  perturbation  surface.  These  waves  can  be  plotted 
in  terms  of  the  flow  variables:  displacement,  strain,  and  particle 
velocity.  If  wave  trapping  is  to  occur,  the  strain  must  tend  to  remain 
localized  at  the  perturbation  surface,  and  only  diffuse  slowly  away. 

That  is,  the  perturbation  acts  as  a  kind  of  piston;  we  are  interested  in 
the  case  when  the  imposed  strain  remains  trapped  at  the  piston. 

This  situation  is  diagrammed  in  Figure  27  In  27(a)  lines  of  con¬ 
stant  strain  are  plotted  in  t  space,  and  in  27(b)  the  corresponding 
profiles  in  h2,  y ^  space  are  shown  at  successive  times.  In  the  situation 
shown,  the  peak  strain  behind  the  perturbation  surface  propagates  more 
slowly  into  the  initially  unstrained  material  than  do  the  lower  strains, 
in  a  diffusion-like  process. 

This  situation  can  also  be  represented  in  terms  of  the  displacement 

field,  u^,  as  shown  in  Figure  28,  which  plots  lines  of  constant  u^ 

radiating  away  from  the  perturbation  surface.  The  dl^/dt  slopes  of  the 

lines  are  the  values  of  C  .  The  gradient  of  the  lines  in  the  h„  direc- 

u^  &  2 

tion  is  -  y  ,  and  the  gradient  in  the  t  direction  is  v^ .  We  see  that 


(8Cu1/8h2>t1 


(9c  /at),  <  o 

U1  h2  ~ 

We  note  that  in  the  region  of  constant  strain,  the  lines  of  constant  u^ 
are  parallel  and  travel  more  slowly  than  the  lines  of  lower  displacement. 
This  corresponds  to  the  diffusion-like  propagation  of  strain  shown  in 
Figure  27(a) ,  from  which  we  see  that 


ay12/9h2  <_  0 


3Y12/,9t  -  0 
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t 


(b)  PROFILES  AT  SUCCESSIVE  TIMES 

MA-64 08-21 

FIGURE  27  EXAMPLE  OF  CASE  WHERE  PLASTIC  WAVE 
TRAPPING  OCCURS 
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NOTE:  Each  path  line  represents  a  constant 
increment  in  u1  from  that  on  the 

ad>aCent  line'  MA-6408-22 

FIGURE  28  PATH  LINES  OF  CONSTANT  DISPLACEMENT  u, 

FOR  CONDITION  OF  LOCALIZATION 


We  next  combine  the  inequalities  (15)  and  (16)  with  the  definitions 
(5)  to  obtain  a  criterion  for  wave  trapping. 

First,  from  the  first  of  the  inequalities  (15),  and  from  (7),  we 
obtain 

y9v/3h  >  v3y/3h  ,  (17) 

where  the  subscripts  have  been  dropped  for  convenience.  Similarly,  from 
the  second  of  the  inequalities  (15)  we  obtain 

y3v/3t<_  v3y/3t  (18) 

Since  y  is  positive,  (17)  and  (18)  can  be  written  as 


-  3y/3  t  >  C  3y/3h 
—  u 


3v/3t  <^Cu(-3v/3h) 


respectively,  where  we  have  used  (7)  to  define  Cu>  Finally,  since  (16) 
shows  that  3y/3  h  is  negative  and  3y/3t  is  positive,  (19)  and  (20)  can 
be  written  as 


C  <  C  , 
y  -  u 


c  <  c  , 

v  —  u 


respectively,  where  we  have  used  (8)  and  (9). 
If  we  now  multiply  (22)  by  C^,  we  obtain 


C  C  >  C 
y  u  —  p 


where  we  have  used  (13)  and  the  fact  that  is  positive.  Relation  (23) 
is  the  one  we  shall  use  as  the  necessary  condition  for  adiabatic  shear 
band  nucleation.  61 


Another  type  of  dynamic  localization  occurs  not  at  the  "piston," 
but  in  a  shock  front.  If  the  strain  increment  introduced  at  the  pertur¬ 
bation  had  propagated  away  as  a  shock,  or  discontinuous  jump  in  strain, 
followed  by  a  decaying  tail,  it  can  be  shown  that  the  inequalities  (15) 
are  unchanged,  but  the  sense  of  inequalities  (16)  are  reversed.  This 
would  cause  a  reversal  of  the  inequality  (23). 

Thus,  if  an  equal  sign  replaces  the  inequality  sign  in  (23),  the 
resultant  equation  defines  a  critical  load  curve,  as  illustrated  in 
Figure  29.  Above  the  curve,  wave  trapping  at  the  perburbation  can 
occur.  Below  the  curve,  perturbations  are  carried  away  from  the  pertur¬ 
bation  to  localize  in  a  shock  font.  That  is,  above  the  curve  deformations 
tend  to  accumulate  at  the  piston  (perturbation),  and  below  the  curve 
they  tend  to  accumulate  in  a  shock  front.  On  the  curve,  the  waves  are 
neither  dispersive  or  convergent,  and  the  perturbation  is  propagated 
unchanged  as  in  an  elastic  material. 

Application  of  the  Criterion  to  Steel 

To  apply  the  above  concept  to  steels,  we  use  the  common  expression 
for  the  work  hardening  curve. 


T 


(24) 


where  T  and  y  are  now  regarded  as  scalar  quantities  expressed  in  terms 

of  suitable  invariants  of  the  stress  and  strain  tensors  and  x  and  n 

o 

are  material  parameters  fitted  to  experimental  data.  We  thus  assume 

1 

isotropic  materials.  The  value  of  C  ,  or  (H/p)2,  is  then  taken  to 


C 

P 


n-1 

Y  2 


(25) 
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29  CRITICAL  LOAD  CURVE  FOR  DYNAMIC  LOCALIZATION 


We  next  write  (23)  as: 


(C  C  ); 
Y  u 


v3  v/3h 

i  >  lm°V  Ui1 

y3y/3h 

-  /  lY2 

(26) 


where  we  have  used  absolute  values  to  avoid  confusion  about  the 

direction  of  the  inequalities.  Table  2  lists  representative  values  of 

* 

T  ,  n,  and  p  for  three  steels,  and  Figure  30  shows  the  corresponding 
o 

critical  load  curves. 


Table  2 


MATERIAL  PROPERTIES  FOR  SEVERAL  STEELS 


Material 


n 


pQ(g/cc) 


Armco  Iron 

5.3  x  109 

0.26 

7.8 

1.3 

X 

10' 

HF-1  steel 

1.57  x  1010 

0.10 

7.8 

1.4 

X 

10* 

4340  steel 

t  o  in10 

3.3  x  10 

0.15 

7.8 

2.5 

X 

io‘ 

RHA 

(1.5"  plate) 

9 

1.2  x  10 

0.08 

7.8 

3.5 

X 

10 

The  data  for  4340,  Rc53  steel  were  generated  at  SRI;  the  RHA  steel 
data  for  the  1.5"  plate  were  taken  from  R.  F.  Benck,  BRL  Memorandum 
Report  No.  2703,  November  1976;  and  the  HF-1  and  Armco  iron  data  were 
estimated  from  data  given  for  similar  steels  in  G.  E.  Dieter, 
Mechanical  Metallurgy,  McGraw-Hill,  N.Y.  (1961),  p.  248. 
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These  critical  load  curves  were  calculated  from  quasistatically 

measured  values  of  T  and  n  in  Eq.  (24).  The  dynamic  work  hardening 

°  4-1 

curves,  corresponding  to  the  strain  rates  of  10  s  often  encountered 

in  shear  banding  experiments,  may  be  significantly  different  for  some 

materials.  Unfortunately,  data  for  these  steels  in  the  high  strain 

13 

rate  region  are  sparse,  although  impact  data  by  Jones  et.  al.  suggest 
that  the  yield  strength  of  4340  steel  is  not  very  rate-dependent.  Hence, 
the  critical  load  curves  plotted  in  Figure  30  may  be  in  error. 


Comparison  With  Experimental  Data 


We  checked  the  plausibility  of  the  above  initiation  criterion  by 
comparing  the  critical  load  curve  for  4340,  R^53  steel  with  calculated 
(C^C^)2  versus  Y  paths  taken  by  material  elements  in  exploding  cylinder 
experiments  that  bracketed  the  threshold  conditions  for  shear  banding. 
Shots  13,  15,  and  16  on  4340,  R^52  steel  all  spanned  the  range  from  no 
shear  banding  in  some  orientation  modes  to  shear  banding  in  other 
orientation  modes. 


In  Shot  13,  material  on  the  inside  surface  of  the  cylinder  exhibited 
incipient  shear  band  striations  in  the  mode  1  orientation,  but  none  in 
the  mode  2  and  3  orientations.  For  this  material  element  we  computed 

i  5 

(C-yC^)  and  y  histories  for  each  oriertation,  and  plotted  the  (C^C^) 
versus  y  paths  in  Figure  30. 


In  contrast  to  Shot  13,  which  showed  only  incipient  mode  1  banding, 
Shot  15  showed  significant  banding  in  the  mode  1  orientation.  A  few 
mode  2  bands  were  found;  these  showed  very  little  growth.  No  mode  3 
loads  were  seen.  The  corresponding  load  paths  are  shown  in  Figure  31. 
Comparison  with  Figure  30  shows  that  in  Shot  15  the  load  paths  for  all 
modes  were  shifted  toward  higher  strains,  but  are  otherwise  very  similar 
to  the  paths  for  Shot  13. 
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FIGURE  31  DYNAMIC  LOAD  PATHS  FOR  SHOTS  15  AND  16 
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Finally,  Shot  16,  performed  at  somewhat  lower  strain  rates,  showed 
two  long  axial  cracks  completely  through  the  cylinder.  One  of  these 
cracks  was  definitely  a  mode  1  shear  band.  We  conclude  that  mode  1 
shear  banding  occurred  in  this  experiment;  the  mode  1  load  path  is 
shown  as  the  dashed  curve  in  Figure  31. 

These  results  are  summarized  in  Figure  32,  which  compares  all  the 
load  paths  for  shear  banded  and  non-shear  banded  material  from  Shots 
13,  15,  and  16. 

Figure  32  shows  that  the  load  paths  for  material  that  nucleated 
shear  bands  are  consistently  above  those  for  which  nucleation  was  not 
observed,  as  expected.  However,  all  the  load  paths  have  portions  above 
the  critical  load  curve,  suggesting  that  nucleation  should  have  occurred 
in  every  case. 

There  are  two  effects  that  could  help  resolve  this  discrepancy. 

First,  as  discussed  above,  the  dynamic  critical  load  curve  may  lie 

above  the  static  curve  shown  in  Figure  32.  Second,  shear  band  growth 

as  well  as  nucleation  is  required  before  the  bands  can  be  observed. 

Our  preliminary  evaluation  of  growth  velocities  suggests  that  the  growth 

law  is  approximately  given  by  Eq.  (2)  in  Section  2,  with  C  being  on 

§ 

the  order  of  10. 

Thus, 

pc 

R/Rq  =  exp  (lOe  )  ,  (27) 

and  the  relative  growth  increases  from  about  3  for  strains  of  0.1  to 
more  than  7  for  strains  of  0.2.  At  still  higher  strains  the  growth 
increases  dramatically.  Thus,  the  observed  transition  from  no  shear 
banding  to  significant  shear  banding  also  corresponds  roughly  to  the 
strain  levels  at  which  significant  growth  begins.  This  is  consistent 
with  the  observed  small  sizes  of  the  mode  2  bands  in  Shot  15. 
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FIGURE  32  COMPARISON  OF  LOAD  PATHS  FOR  SHEAR  BANDED  AND 


NON-SHEAR  BANDED  MATERIAL 
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Discussion  of  Postulated  Initiation  Criterion 

The  postulated  initiation  criterion  based  on  plastic  wave  trapping 
appears  at  this  stage  to  be  plausible,  but  requires  further  investigation. 
In  future  work  we  will  attempt  to  obtain  dynamic  stress-strain  data  for 
construction  of  dynamic  critical  load  curves.  We  will  also  attempt  to 
further  improve  the  growth  law  by  analysing  additional  shear  band  size 
data. 


70 


6  SIMPLIFIED  TESTS  TO  OBTAIN  SHEAR  BANDING  DATA 


We  examined  a  number  of  experimental  techniques  that  have  been 
used  or  proposed  for  use  in  studying  shear  banding  phenomena.  We 
evaluated  the  techniques  for  their  ability  to  provide  data,  at  minimal 
time  and  expense,  from  which  we  could  obtain  the  parameters  required  as 
input  to  our  shear  band  model.  These  parameters  are  basically  a  thres¬ 
hold  or  initiation  parameter  (to  specify  if  the  loading  conditions  allow 
shear  bands  to  form),  a  nucleation  parameter  (to  specifiy  the  number  or 
density  of  bands  that  have  been  formed),  and  a  growth  parameter  (to 
specify  the  growth  velocity  of  the  bands).  These  three  parameters  are 
functions  of  loading  conditions  (e.g.,  the  plastic  shear  strain  rate), 
and  the  amount  of  damage  is  related  to  the  duration  of  this  applied  load. 
To  obtain  the  shear  hand  parameters,  one  needs  an  experimental  technique 
in  which  magnitude  and  duration  of  the  loading  can  be  varied  more  or 
less  independently. 

In  addition  to  the  contained  fragmenting  cylinder  technique  described 
in  detail  in  Section  3,  there  are  a  number  of  other  experimental  tech¬ 
niques  for  producing  shear  bands  in  the  laboratory.  They  can  be  divided 
into  three  broad  classes  of  techniques:  dynamic  punch  tests,  dynamic 
torsion  tests,  and  high-speed  machining  tests. 

The  dynamic  punch  tests  include  all  forms  of  impact  or  penetration 

experiments,  in  particular  those  using  a  flat-nosed  cylindrical  projectile 

wherein  shear  banding  in  the  target  parallel  to  the  impact  direction  is 

14 

the  cause  of  target  failure.  The  most  promising  technique  of  this 

type  appears  to  be  the  "stepped"  projectile  technique,  described  by 

15  33 

Wingrove  and  shown  schematically  in  Figure  .  The  magnitude  and 

duration  of  the  load  can  be  varied  independently  by  altering  the  impact 
velocity  and  step  length,  respectively.  This  technique  was  used  success¬ 
fully  to  study  shear  banding  in  aluminum  targets  impacted  by  hardened 
steel  penetrators.  However,  the  shear  load  is  concentrated  in  such  a 


PROJECTILE 


TARGET  SUPPORT 

(a)  (b) 


(O 

MA-6408-26 


FIGURE  33  SCHEMATIC  OF  WINGROVE’S  STEPPED  PROJECTILE  TECHNIQUE  15 

Stepped  projectile  is  fired  at  supported  target  (a).  Step  penetrates  target  (b| 
until  depth  of  penetration  equals  step  length.  Then  projectile  body  impacts 
target,  "freezes"  step  penetration,  and  shears  out  large  plug  (c)  for  soft 
recovery. 
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narrow  region  that  it  is  extremely  difficult  to  measure  or  calculate 
the  stress  or  strain  states  existing  in  that  region.  Therefore, 
although  it  may  be  possible  to  obtain  a  predictive  capability  for  flat¬ 
nosed  penetration,  one  could  not  obtain  the  parameters  necessary  for  a 
general  shear  band  model  that  would  be  useful  in  other  geomtries. 
Furthermore,  this  technique  may  not  be  suitable  for  use  with  an  armor 
steel  as  the  target  material  because  of  large  deformation  in  the  pene- 
trator.  Nevertheless,  this  technique  is  valuable  in  identifying  pheno¬ 
menology  for  armor  penetration  applications. 

Dynamic  torsion  tests  with  thin-walled  tubular  specimens'*^  have 

been  used  to  generate  shear  bands  in  steel  at  shear  strain  rates  of 
2 

«10  /sec.  The  loading  magnitude  and  duration  can  be  varied  by  changing 
the  angular  velocity  and  maximum  angular  stroke  of  the  torsional  testing 
machine.  However,  the  maximum  strain  rates  of  machines  of  this  type 
are  orders  of  magnitude  lower  than  the  rates  found  in  exploding  rounds 
or  penetration  experiments.  This  is  a  serious  drawback,  because  shear 
bands  are  very  sensitive  to  strain  rate. 

High-speed  machining  of  metal  has  been  used  to  study  shear  band 
17 

formation,  and  may  in  the  future  be  used  to  measure  the  temperatures 

18 

generated  during  the  shear  band  process.  However,  the  extreme  load 
concentration  in  the  region  of  band  formation  is  also  the  major  drawback 
to  use  of  this  technique  for  generating  shear  band  model  parameters. 

Our  present  conclusion  is  that  the  contained  fragmenting  cylinder 
tests  provide  the  most  unambiguous  data  for  shear  band  modeling.  The 
loads  and  strains  are  evenly  and  smoothly  distributed  over  the  specimen, 
thus  allowing  use  of  a  maximum  volume  of  material  in  a  single  test. 
Furthermore,  the  simple  geometry  lends  itself  to  computer  modeling. 

For  these  reasons,  we  recommend  that  this  method  be  used  to  obtain  shear 
banding  parameters  for  any  material  that  can  be  obtained  in  tubular  form. 
However,  these  tests  should  be  supplemented  by  tests  that  approximate 
the  geometrical  constraints  of  the  applications  of  interest.  For 
example,  for  armor  penetration  with  long  rod  penetrators,  dynamic  punch 
tests  are  valuable  for  identifying  the  key  phenomenological  processes. 


7  DISCUSSION  AND  RECOMMENDATIONS 


The  experimental  program  has  fulfilled  our  objective  of  providing 
a  data  base  for  4340  steel  that  allows  improvement  of  the  SNAG  model. 
Data  from  expanding  cylinder  tests  now  exist  for  three  hardnesses  of 
4340  steel,  spanning  shear  banding  levels  of  damage  from  incipient 
shear  banding  to  complete  fragmentation.  For  one  hardness  (Rc40) , 
quantitative  number,  size,  and  orientation  distributions  for  the  shear 
bands  were  measured  and  correlated  with  calculated  strain  and  stress 
histories.  These  correlations  suggest  that  strain  rate  is  indeed  the 
key  loading  parameter  that  governs  the  number  and  growth  of  the  shear 
bands.  Furthermore,  our  derived  nucleation  threshold  condition  suggests 
that  the  shape  of  the  dynamic  work  hardening  curve  is  crucial  for 
determining  the  strain  and  strain  rate  conditions  required  to  nucleate 
the  bands. 

The  above  results  indicate  that  the  following  future  work  is  highly 
desirable: 

•  Complete  the  quantitative  data  analysis  by  measuring  the 
shear  band  size  and  orientation  distributions  for  the 
Rc52  cylinders. 

•  Measure  the  static  and  dynamic  stress-strain  curves  for 
all  three  hardnesses.  According  to  Jones  et.  al.,  (Ref.  13) 

4340  steel  may  not  be  excessively  rate-dependent,  but  this 
should  be  checked,  at  least  at  rates  attainable  on  standard 
tensile  testing  machines. 

•  Continue  to  analyze  the  data  with  iterative  computations 
to  obtain  the  SNAG  parameters  as  functions  of  the  strain 
hardening  curve. 

•  Perform  dynamic  punch  tests  and  computer  simulations  for 
targets  with  each  of  the  three  hardness  levels  to  see  if 
the  SNAG  parameters  developed  from  expanding  cylinder 
tests  will  correctly  predict  and  rank  the  phenomenology 
of  shear  banding  in  geometries  applicable  to  armor  pene¬ 
tration. 


The  above  tasks  are  needed  to  fully  use  the  experimental  data 
base  generated  in  this  program;  they  will  provide  an  understanding  of 
how  the  shear  banding  process  depends  on  key  material  properties  for 
one  material.  Such  an  understanding  is  required  to  allow  extrapolation 
to  other  materials  such  as  those  used  for  armor  and  penetrators. 
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1  ATTN:  Technical  Library 
1  Dr.  E.  W.  Ross 

1  ORONA-UE ,  Dr.  L.  A.  McClaine 

1  Mr.  .J.  Flanagan 

f<»*i"ander ,  U.  S.  Army  Tank 'Automotive  Research  and  Development 
Cornnand,  Warren,  Michigan  48090 
1  ATTN:  DRDTA-UL ,  Technical  Library 

1  Mr .  V .  Paqano 

Cormande**,  U.  s.  Army  Armament  Research  and  Development 
Command,  Dover,  New  Jersey  07801 
1  A1TN:  Technical  Library 
1  DP0AP-5CM,  J.  D.  Corrie 

1  Dr.  J.  Brinkman 

1  Or.  J.  Miller 

1  Mr.  J.  Beetle 

1  Dr.  J.  F raiser 

1  Mr.  Harry  E.  Pebly,  Jr.,  PLASTEC,  Director 

1  Mr.  T.  Knight 

1  Mr.  .1.  Mershkowitz 

1  Mr.  R.  Knnrad 

1  Mr.  W.  Pryar 

1  Mr.  J.  Pearson 

1  G.  Randers-P*hrson 

1  Mr  A  Garcia 

Commander,  Rock  Island  Arsenal ,  Rock  Island,  Illinois  61299 

1  4HN;  Mr.  R.  Nelson 
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Director,  U.  S.  Army  Ballistic  Research  Laboratory, 

Aberdeen  Proving  Ground,  Maryland  21005 
1  ATTNi  Dr.  A.  Dietrich 

1  Dr,  J.  Suckling 

1  Dr.  J.  Zukas 

1  Dr.  R.  Vi  tali 

1  Or  G.  L.  F i 1  bey 

1  Dr.  W.  Gillich 

1  Dr.  G.  Moss 

1  OROAR-TSB-S  (STINFO) 

Commander  Harry  Diamond  Laboratories,  2800  Powder  Mill  Road, 
Adelphi,  Maryland  20783 
1  ATTN:  Technical  Information  Office 

Commander,  Watervliet  Arsenal ,  Watervliet,  New  York  12189 
1  ATTN:  Dr,  T.  Davidson 

1  Mr.  D.  P.  Kendall 

1  Mr  J.  f.  Thrc-T- 

Commander,  U.  S.  Army  Foreign  Science  and  Technology  Center, 
220  7th  Street,  N.  E.,  Charlottesville,  Virginia  22901 
1  ATTN:  Mr.  Marley,  Military  Tech 

Chief,  Benet  Weapons  Laboratory,  LCWSL,  USA  ARRADCOM, 
Watervliet  Arsenal,  Watervliet.  New  York  12189 
1  ATTN;  DRDAR-LCB-TL 

Director,  Eustis  Directorate,  U.  S.  Army  Air  Mobility 
Research  and  Development  Laboratory,  Fort  Eustis, 

Virginia  23604 

)  ATTN:  Mr.  J.  Robinson,  DAVDl-F-MOS  (AVRflDCOM) 

Commander ,  USACOC  Air  Defense  Aqency,  Fort  Bliss, 

Texas  79916 

1  ATTN:  Technica1  Library 

Cownander,  IJ.  S.  Army  Fnnineer  School,  Fort  Belvoir, 

Virginia  22060 
1  ATTN:  Library 

Conmander,  U.  S.  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  Mississippi  39180 
1  ATTN:  Research  Center  Library 

Comnander,  U.  S.  Army  Mobility  Equipment  Research  and 
Development  Center,  Tort  Belvoir,  Virginia  22060 
1  ATTN:  Mr.  william  H.  Baer,  Materials  Research  Laboratory 

Director,  Structural  Mechanics  Research,  Office  of  Naval 
Research,  800  North  Quincy  Street,  Arlington,  Virginia  22203 
1  ATTN:  Dr.  N.  Perrone 

Naval  Research  Laboratory,  Washington,  D.C.  20375 
1  ATTN:  C.  D.  Rearhem,  Head,  Adv.  Mat’ Is  Tech  Br.  (Code  6310) 

1  Dr.  J.  M.  KrafU  -  Code  8*30 

1  Dr.  P.  P.  Puzak 

1  Dr.  Jim  C.  I .  Chang 

Naval  Ship  Research  and  Development  Laboratory, 

Annapolis,  Maryland  2140? 

1  ATTN:  Dr.  H.  P.  Chu 

Chief  of  Naval  Researrh,  Arlington.  Virginia  22217 
1  ATTN:  Code  471 

Naval  Surface  Weapons  Center,  Oahlgren,  Virginia  2244R 
1  ATTN:  W.  Mannschrec V 
1  W.  Mock 

1  0.  Brunson 

1  Mr.  W.  H.  Holt 

1  Technical  Library 
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Naval  Surface  Weapons  Center,  Silver  Sprinq.  Maryland  20910 
1  ATTN:  Or.  P..  Crowe 
1  Technical  Library 

Naval  Weapons  laboratory,  Washington,  D.C.  20  390 
1  ATTN:  H.  W.  Pomine,  Mail  Stoo  103 

Naval  Weapons  Center,  China  Lake,  California  93555 
1  ATTN:  Mr.  M.  Backman 

1  Dr.  J.  Pearson 

Air  Force  Materials  Laboratory,  Wriqht-Patterson  Air  Force 
Base.  Ohio  45433 

2  ATTN,  AFML  (MXE),  E.  Morrissey 

1  AFML  (LC) 

1  AFML  (LLP),  D.  M.  Fomev,  Jr. 

1  AFML  (LNC),  T,  J.  Reinhart 

1  AFML  fMBC),  Mr.  Stanley  Schulman 

1  AFFOL  (TB),  Dr.  J.  C.  Halpin 

1  Dr.  T.  Nicholas 
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Air  Force  Base,  Ohio  45433 
1  ATTN-  AFFOL  (fbs),  C.  Wallace 
1  L i bra  r  v 

National  Aeronautics  and  Space  Administration, 

Washington,  D.C.  20546 
1  ATTN:  Mr.  B.  G.  Achhammer 
1  Mr.  G.  C.  Deutsch  -  Code  RW 

National  Aeronautics  and  Space  Administration,  Marshall 
Space  Fliqht  Center,  Huntsville,  Alabama  35812 
1  ATTN:  R.  J.  Schwinqhamer,  EH01 ,  Dir.,  MSP  Lab 
1  Mr.  W.  A.  Wilson,  EH41 ,  Bldg.  4612 

National  Bureau  of  Standards,  U.  S.  Department  of  Cornnerce, 
Washington.  D.C.  20234 
1  ATTN-  J.  A.  Bennett 

T  Mechanical  Properties  Data  Center,  Belfour  Stulen  Inc., 
13917  W.  Ray  Shore  Drive,  Traverse  City,  Michiqan  49684 

Honevwell,  I  nr.,  Pesearch  Dept.,  SSRC,  Mail  Stop  A3  340* 
2345  Walnut  Strpet,  Saint  Paul,  Minnesota  55113 
1  ATTN-  Mr.  Gene  Fisher 
1  Or.  Gordon  Johnson 

1  Prof.  I.  Fyfe,  University  of  Washington,  Department  of 
Aeronautics  and  Astronautics,  Seattle,  Washington  98105 

Stanford  Research  Institute,  Menlo  Park,  California  94025 
1  ATTN:  Dr.  D.  Curran 

1  Dr.  L.  Seaman 

1  Dr.  D.  Snockey 

1  Dr.  G.  Abrahamson 

1  Dr.  0.  Erlich 

1  Dr.  R.  CaTigiuri 

Systems,  Science  and  Software,  P.  0.  Box  1620,  La  Jolla, 
California  9?D*7 
1  ATTN  Dr.  P.  ^edowick 

Sandu  Corpora on,  Mrtland  East,  Albuquerque, 

New  Mexico  D/115 
1  ATTN  [)r.  B.  Butrhpr 

1  Dr.  A.  r>tevpns 

1  hr.  w.  Uprrtnann 

1  Dr.  L.  Prrtholf 

1  Pr .  5 .  Rpn^lpy 

Southwest  Research  Institute,  8500  Culebra  Road, 

Sar  Antonio,  Texas  7P?84 
1  AT'V  Prnf.  !j,  Lindholm 
1  Dr.  W.  T.akPr 

1  Dr.  Geo rap  0  Irwin.  Department  of  Mechanical  Engineering, 
unsprfts  nf  v^rvl *nd ,  College  Park,  Maryland  20742 


No.  of 

Copies  To 

Battelle  Columbus  Laboratories ,  50^  king  Avenue,  Colunfcus, 

Ohio  43201 
1  ATTN:  E.  Hulbprt 

Brown  University,  Providence,  Rhode  Island  02912 
1  ATTN:  Prof.  J.  R.  Rice 

1  Prof.  W.  N.  Findley,  Division  of  Engineering,  Box  D 

1  Prof.  R.  Clifton 

1  Prof.  J.  Duffey 

1  Carnegie-Mel Ion  University,  Department  of  Mechanical 

Engineering,  Schenley  Park,  Pittsburgh,  Pennsylvania  15213 
1  ATTN:  Dr.  J.  L.  Swedlow 

1  Prof.  J.  D.  lubahn,  Colorado  School  of  Mines,  Golden, 

Colorado  80401 

1  Prof.  J.  Dvorak,  Civil  Engineering  Department, 

University  of  Utah,  Galt  Lake  City,  Utah  84112 

George  Washington  University,  School  of  Engineering  and 
Applied  Sciences,  Washington,  D.C.  20052 
1  ATTN:  Dr.  H.  Liebowitz 

Terra  Tek,  Univprsity  Research  Park,  420  Wakara  Way, 

Salt  Lake  City,  Utah  84108 
1  ATTN:  Mr.  S.  Green 
1  Dr.  A.  Jones 

Drexel  University,  Department  of  Materials  Science, 

Market  and  Chestnut  Streets,  Philadelphia,  Pennsylvania  19104 
1  ATTN:  Prof.  H.  C.  Rogers 

Washington  State  University,  Department  of  Physics, 

Pullman,  Washington  99163 
1  ATTN:  Prof.  G.  Duvall 
1  Prof.  R.  Fowles 

Massachusetts  Institute  of  Technology,  Cambridge, 

Massachusetts  02139 

1  ATTN:  Prof.  F .  A.  MrClintOck,  Pnom  1-304 

1  Prof.  R.  M.  Pelloux 

1  Prof.  T.  h.  w.  Pi  an,  Department  of  Aeronautics 

and  Astronautics 

1  Prof.  A.  S.  Argon,  Room  1-306 

1  Dr.  Greqory  Olson,  Dept,  of  Metallurgy 

1  Dr.  0.  Parks,  Dent,  o'  Mech.  Eng. 

Syracuse  University.  Department  of  Chemical  Engineering 
and  Metalluroy,  409  link  Hall,  Syracuse,  New  York  13210 
1  ATTN:  Dr.  V.  Weiss,  Meta  11 urqi ca l  Research  Labs.,  Bldg.  D-6 

1  Prof.  W.  Goldsmith,  Department  of  Mechanical  Engineering. 
University  of  California.  Berkeley,  California  94720 

University  of  California,  Los  Alamos  Scientific  Laboratory, 
Los  Alamos,  New  Mexico  87544 
1  ATTN:  Dr.  R.  Karcp 
1  Dr.  J.  M.  Waist 

1  Technical  library 

1  Prof.  A.  J.  McEvily,  Metallurgy  Department  U- 1 36 *  University 
of  Connecticut.  Storrs,  Connecticut  06268 

1  Prof.  D.  Drucker.  Dean  of  School  of  Engineering,  University 
of  Illinois,  Champaiqn,  Illinois  61820 

University  of  Illinois,  Urbana,  Illinois  61801 
1  ATTN:  Prof.  T.  ■).  Dolan,  Department  of  Theoretical  and 
Appl ied  Mechanics 

1  Prof.  T.  la*dnp,-.  Dept.  cf  Theoretical  and  Applied 

Mechanics 

1  Dr.  M.  L.  Williams,  Doan  of  Eoqinepring,  240  Benedum  Hall. 
University  of  Pi tt sbw» tih ,  Pittsburg,  Pennsylvania  1526D 

1  Prof.  A.  k.obayash:,  Department  of  Mechanical  Engineering, 
FU-10,  University  of  Washington.  Seattle,  Washington  98195 
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1  R.  Sedgwick,  Systems  Science  Softwave,  La  Jolla, 

California  92037 

1  Prof.  W.  f.  Riley,  Dept,  of  Eng.  Mech.,  Iowa  State 
University,  Ames,  Iowa  50010 

1  Prof.  I.  Fyfe,  Dept,  of  Ae**o.  Eng.,  University  of 
Washington,  Seattle,  Washington  98105 

State  University  of  New  York  at  Stony  Brook,  Stony  Brook, 
New  York  11790 

1  ATTN:  Prof.  Fu-Pen  Chianq,  Department  of  Mechanics 

Lawrence  Livermore  Laboratory,  Livermore,  California  94550 
1  ATTN:  Or.  R.  To! and 

1  Dr.  M.  Wilkins 

1  Dr.  M.  Guinan 
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Denver  Research  Institute,  2390  South  University  Boulevard. 
Oenver,  Colorado  80210 
1  ATTN:  Dr.  «  Recht 

1  Dr.  C.  Hoggart 

Director,  Army  Materials  and  Mechanics  Research  Center, 
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2  ATTN:  DRXMR-PL 

1  DRXMR-PR 

1  DRXMR-PD 

1  DRXMR-AP 

26  DRXMR-TE,  Mr.  J.  Mescal  1 
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Prelimmary  computer  simulations  were  performed  to  determine  plastic  shear  strain  rate  Preliminary  computer  simulations  were  performed  to  determine  plastic  shear  strain  rate 

histories  in  the  direction  of  the  predominant  sheer  banding  orientation  mode,  and  the  com-  histories  in  the  direction  of  the  predominant  shear  banding  orientation  mode,  and  the  com¬ 
putations  were  compared  with  experimental  data  to  yield  first  estimates  for  the  values  of  the  putations  were  compared  with  experimental  data  to  yield  first  estimates  for  the  values  of  the 


